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ABSTRACT 

Context. Current knowledge suggests that the dust-driven wind scenario provides a realistic framework for understanding mass loss 
from C-rich AGB stars. For M-type objects, however, recent detailed models demonstrate that radiation pressure on silicate grains is 
not sufficient to drive the observed winds, contrary to previous expectations. 

Aims. In this paper, we suggest an alternative mechanism for the mass-loss of M-type AGB stars, involving the formation of both 
carbon and silicate grains due to non-equilibrium effects, and we study the viability of this scenario. 

Methods. We model the dynamical atmospheres and winds of AGB stars by solving the coupled system of frequency-dependent 
radiation hydrodynamics and time-dependent dust formation, using a parameterized description of non-equilibrium effects in the gas 
phase. This approach allows us to assess under which circumstances it is possible to drive winds with small amounts of carbon dust 
and to get silicate grains forming in these outflows at the same time. 

Results. The properties of the resulting wind models, such as mass loss rates and outflow velocities, are well within the observed 
limits for M-type AGB stars. Furthermore, according to our results, it is quite unlikely that significant amounts of silicate grains 
will condense in a wind driven by a force totally unrelated to dust formation, as the conditions in the upper atmosphere and wind 
acceleration region put strong constraints on grain growth. 

Conclusions. The proposed scenario provides a natural explanation for the observed similarities in wind properties of M-type and 
C-type AGB stars and implies a smooth transition for stars with increasing carbon abundance, from solar-composition to C-rich AGB 
stars, possibly solving the long-standing problem of the driving mechanism for stars with C/O close to one. 
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While alternative wind mechanisms have been discussed ear- 
lier, observations indicate that winds of O-rich AGB stars be- 



^ 1. Introduction these models show that the Fe-content of the grains has to be 

. extremely low, as the radiative equilibrium temperature of sil- 

Low- and intermediate-mass stars are known to lose a significant icate infj increasef , j ^ increasi inclusion of Fe> 

g . fraction of their mass through slow massive winds during the which leads tQ eyen ^ in ^ than iousl as _ 
Asymptotic Giant Branch (AGB) phase of their evolution. The sumed the ingre dient that makes C-rich mod- 
basic scenario for this process is that of pulsation-enhanced dust- dg mQre successful than evel . j e frequency-dependent radiative 
driven winds: shock waves created by stellar pulsation lead to a ntfj dufjt driyen windfj ^ starg wkh < 1. 
dense, cool extended stellar atmosphere, allowing for efficient 
dust formation. The grains are accelerated away from the star by 
radiation pressure, dragging gas along. The composition of the 
grains is determined by the elemental abundances in the stellar have ver y Slmllar t0 their C - nch counterparts. Furthermore, the 
atmosphere, with C-rich stars forming carbon grains, and O-rich ver y fact that Slllcate g rains m formed ln such outflows P uts 
environments producing silicate particles. rather stnct constraints on the possible mechanisms, as will be 
This scenario seems to work well for the case of C-rich AGB discussed in Sec. ED This leads us to propose a scenario, where 
stars, as demonstrated by comparison of detailed self-consistent Wlnds of (and probably S-) type AGB stars are actually 
dynamical models with observations as diverse as low-resolution dnven b ? sma11 amounts of carbon g rains ' wlth Slllcates form - 
IR spectra (e.g., Hofner et al. 2003, Gautschy-Loidl et al. 2004) ln 8 as a by-product. This requires that a certain fraction of the 
and line profile variations of CO vibration-rotation lines (e.g., carbon atoms in the dust formation region are available for grain 
Nowotny et al. 2005). formation, i.e. not bound up in CO. 

Recent work on similar models for stars with C/0< 1, This scenario implies deviations from chemical equilibrium 

however, faces serious problems (Woitke 2006; Hofner 2007). (CE) in the gas phase, probably brought about by strong atmo- 

Frequency-dependent dynamical models demonstrate that the spheric shock waves, which seems not too far-fetched in the 

opacities of silicate grains are too low to drive a wind, even light of recent observations showing the existence of molecules 

when assuming rather favorable grain compositions. In addition, typically expected in O-rich environments in C-rich stars (e.g., 

Schoier et al. 2006), and vice versa. Current models of non- 

Send offprint requests to: S. Hofner equilibrium gas-phase chemistry by Cherchneff (2006) predict 
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rather small deviations from CE (less than what we assume be- 
low), but these models are computed for significantly higher 
densities (which should favor more CE-like compositions), and 
without taking condensation into account. As an inclusion of de- 
tailed non-CE in the gas phase in our wind models is well be- 
yond the scope of this first investigation, we use the abundance 
of available carbon as a parameter of the models. This approach 
allows us to test under which circumstances it is possible to drive 
winds and to get a significant amount of silicate grains forming 
in the outflows. 



2. The hydrodynamical model 

Our hydrodynamical models predict mass-loss rates and wind 
velocities of AGB stars, as well as the amount of dust formed 
in the winds, by treating in detail the atmosphere and the cir- 
cumstellar environment around pulsating long-period variable 
stars. This is done by solving the coupled system of frequency- 
dependent radiation hydrodynamics and time-dependent dust 
formation (cf. Hofner et al. 2003). The calculations presented 
here are based on opacity sampling data of molecular opacities 
at 64 frequency points between 0.25 and 25 fim. 

In contrast to our earlier models, we have implemented a 
time-dependent description of grain growth for silicate grains 
(Hofner et al., in prep.), and we focus on the formation of pure 
forsterite particles (Mg2SiC>4), as non-grey effects will force a 
low Fe-content (see Sect. Q]). While treating the growth of grains 
in full non-equilibrium, we do not consider the nucleation of new 
grains from the gas phase but assume the presence of seed nuclei 
(a given number per H atom) at the point where grain growth 
becomes possible. Opacity data for the relevant grain materials 
is taken from Jager et al. (2003; forsterite), Rouleau & Martin 
(1991; amorphous carbon sample AC2) and Jager et al. (1998; 
amorphous carbon sample cel400). 

All models discussed here have solar elemental abundances 
(i.e. C/O « 0.5), but we assume that a certain fraction of car- 
bon is not bound in CO in the region relevant for grain growth. 
This fraction is parameterized by a factor Fc = £c/sc, where sc 
is the abundance of carbon available for dust formation. This is 
admittedly a very simple-minded treatment of non-equilibrium 
in the gas phase, but it allows us to assess, how much free car- 
bon is needed to drive a wind for given stellar parameters. Note 
also, that only a fraction of this fraction of all carbon atoms will 
actually condense into grains (see Tab. Q]and discussion below). 

3. Results and discussion 

We have calculated a number of models with different combina- 
tions of stellar parameters, fractions of free carbon (.Fc), and 
pulsation amplitudes, using two different sets of opacity data 
for amorphous carborQ. The model parameters and the result- 
ing wind properties are summarized in Tab. [1] As expected, the 
fraction of free carbon required to obtain a dust-driven wind de- 
creases with increasing luminosity and decreasing effective tem- 
perature. Also, in accordance with previous results on C-rich 
wind models, the mean degree of condensation of free carbon 
fc is never close to one in these models, and the degree of con- 
densation of all carbon (obtained by multiplying fc with Fc) is 
rather low in certain cases. 



1 This term actually covers a variety of materials with different mi- 
croscopic structures and different optical properties (see discussion in 
Andersen et al. 2003). 



What may seem more surprising at a first glance, is the fact 
that the degree of condensation for Si (and, correspondingly, the 
amount of silicate grains formed in the wind) is quite low in most 
models, and actually decreases with conditions that make the 
formation of carbon grains more efficient (e.g., higher pulsation 
amplitude, or a larger fraction of free carbon, for given stellar pa- 
rameters). The low silicate dust abundance is not a matter of tem- 
peratures being too high to allow for condensation, but rather a 
matter of grain growth timescales: as soon as enough carbon dust 
is formed to accelerate the gas away from the star, the condensa- 
tion of dust (and in particular of silicate grains, handicapped by 
lower abundances of critical elements like Si and Mg) becomes 
a race against falling densities in the outflow, a phenomenon that 
finally stops condensation, both for silicates and carbon grains 
(leading to fc < 1), despite falling gas and grain temperatures. 

Therefore it is actually quite remarkable that silicate grains 
can form at all in a wind driven by some other force (e.g. radia- 
tion pressure on another, unrelated dust species as in this case), 
as the rapid slowing down of grain growth in an outflow puts 
tight constraints on the wind mechanism. If the silicate grains 
do not form simultaneously with (or maybe even prior to) the 
driving dust species (carbon), i.e. at the same distance from (or 
maybe even closer to) the star, at rather high gas densities, there 
is virtually no chance of their forming at all. 

In order to understand how grains consisting of such dif- 
ferent materials as amorphous carbon and silicates may form 
at basically the same distance from the star, we have to take a 
look at two types of temperatures: the condensation tempera- 
ture of a material (setting the threshold for grain formation) and 
the grain temperature at a given distance from the star, which, 
in practice, will be equal to the radiative equilibrium tempera- 
ture corresponding to the grain material. At the densities pre- 
vailing in the pulsating atmosphere, amorphous carbon grains 
will typically start to form around 1 500-1 800 K while the forma- 
tion of forsterite requires temperatures as low as 1000-1 200 K. 
However, due to their optical properties (opacity decreasing with 
wavelength around stellar flux maximum), carbon grains will be 
hotter than the black body temperature at a given distance (cf. 
Hofner et al. 2003, Fig. 6), while forsterite grains will be cooler 
(opacity increasing with wavelength in the critical regime)0This 
means that, at the same distance from the star, forsterite grains 
will be significantly cooler than amorphous carbon grains, and 
the critical condensation temperature is reached approximately 
in the same region for both types of dust (Fig.[TJ. 

The strong sensitivity of a non-driving dust species (i.e. sil- 
icate grains in the present case) to the conditions in the wind 
acceleration region is illustrated by a comparison of models cal- 
culated with two different opacities for the carbon grains, but 
otherwise identical parameters. The data of Jager et al. (1998; 
cel400) has a steeper slope in the region around 1 /vm, than many 
other types of amorphous carbon (cf. Fig. 2 in Andersen et al. 
2003) leading to a higher grain temperature at a given distance. 
Therefore, the formation of carbon grains (and, consequently the 
wind acceleration region) is shifted slightly outwards in these 
models compared to the ones based on the Rouleau & Martin 
(1991) data, improving the conditions for the formation of sili- 
cate grains. This results in an increased degree of condensation 



2 At the relevant wavelengths, the opacity of amorphous carbon 
grains is roughly proportional to \/A. Therefore the carbon grains ab- 
sorb and emit radiation preferentially at short wavelengths. This prop- 
erty in combination with the radiative equilibrium condition requires 
the grains to be hotter than the radiation temperature in order to emit as 
much energy as they absorb. 
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Table 1. Model parameters and input data: luminosity L*, effective temperature (initial model; mass M* = 1M for all models); 
fraction of free carbon Fq, period P and velocity amplitude of piston Au p (luminosity factor fi = 2.0 for all models, see Gautschy- 
Loidl et al. 2004); opacity data for amorphous carbon kq'- RM = data from Rouleau & Martin (1991, sample AC2), J4 = data from 
Jager et al. (1998, sample cel400). Resulting wind properties: mass loss rate M, mean velocity at the outer boundary (u), degree of 
condensation of silicon fsu of free carbon fc and of all carbon f„ ot . 
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of Si, while the mass loss rate and wind velocity (governed by 
the carbon grains) are basically unchanged (cf. Tab. [T]0. These 
results indicate that it may be very difficult to form silicate grains 
in a wind driven by a force totally unrelated to dust formation, 
as the conditions in the upper atmosphere and wind acceleration 
region put strong constraints on grain growth. 

In addition to these theoretical considerations, several ex- 
isting observations can be interpreted in favor of the proposed 
scenario: Schoier et al. (2006) have recently demonstrated the 
existence of molecules typical of O-rich environments in C- 
rich AGB stars, a finding that is hard to reconcile with stan- 
dard chemical equilibrium models. Interferometric and spectro- 
scopic observations indicate the existence of very dense layers 
of water vapor in O-rich stars (e.g., Weiner 2004, Ohnaka 2005, 
Justtanont 2005), potentially requiring abundances of H2O in ex- 
cess of what can be expected from chemical equilibrium (where 
an amount of oxygen equal to the total amount of carbon is 
bound in CO). Furthermore, while the wind properties of C-rich 
and O-rich AGB stars are very similar, there is an average ten- 
dency of slightly higher outflow velocities at given mass loss 
rate for C-rich stars (cf., e.g., Olofsson 2004). This could be a 
natural consequence of all these winds being driven by carbon 
grains, but with an increasing amount of carbon being available 
for condensation, going from O-rich to C-rich objects. 

4. Conclusions 

We have investigated the viability of a scenario where winds of 
O-rich AGB stars are driven by small amounts of carbon grains. 
Our approach allows us to estimate both lower and upper lim- 
its for the required non-equilibrium effects in the gas, with one 
end of the spectrum being defined by the minimum amount of 
carbon dust necessary to drive an outflow, and the other by the 
decreasing amounts of silicates being formed in the wind with 
increasing availability of carbon. Our model shows that signif- 
icant amounts of forsterite grains (recently found in comets by 
STARDUST) can form in such winds. 

The proposed scenario provides a natural explanation for the 
observed similarities in wind properties of O-rich and C-rich 
cool giants and implies a smooth transition for stars with in- 
creasing carbon abundance, from M-type to C-type AGB stars, 



probably solving the long-standing problem of the driving mech- 
anism for stars with C/O close to one (S-type AGB stars; cf., e.g., 
Ramstedt et al. 2006). 

Crucial, if non-trivial, tests of the scenario can be provided 
by observations which can assess the relative abundances of 
molecules in the critical dust formation region, in order to check 
for non-equilibrium effects in the gas phase. A direct detection 
of carbon grains in the spectra is rather unlikely: the optical 
properties of amorphous carbon lead to remarkably featureless 
IR spectra. SiC, on the other hand, which is often used as an 
indicator for dust in C-rich stars due to its pronounced feature 
around 1 1/im will hardly be produced in sufficient quantities to 
be detected above the background of amorphous carbon and the 
nearby silicate feature. Considering the relative opacities of SiC 
and silicates, this would require at least as much Si being bound 
in SiC as in silicates, which seems unlikely in an environment 
with C/O < 1. 
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3 Note, again, that the degree of condensation of Si is a direct mea- 
sure of the amount of silicate grains formed in the wind. 
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Fig. 1. Snapshot of the radial structure of the model with L* = 7000 Lq, T+ = 2700 K and opacity data for amorphous carbon by 
Jager et al. (1998; sample cel400). Top row: gas velocity and density; middle row: radiative equilibrium temperature for forsterite 
and carbon grains; bottom row: degree of condensation for Si and C. Radius in units of the stellar radius. 



